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Abstract 

Although the existence of asymmetrical profile of [OIII]A5007 has been discovered 
for ages, its filiation and physics are poorly understood. Two new spectra of I ZWI 
taken on Nov 16, 2001 and on Dec 3, 2002 were compared with the spectra taken 
by BG92. Following results are obtained. l)The certain variations of broad [OIII] 
during about 10 years separating the observations are identified. The inferred length 
scale of broad [OIII] emitting region ranges from 0.3pc to 3pc. By assuming a 
Keplerian motion in emitting region, the material emitting broad [OIII] is likely to 
be located at transient emission line region, between BLR and NLR. 2)We find a 
positive relation between the Fell emission and flux of H/?(or continuum). On the 
other hand, the parameter RFe decreases with ionizing continuum marginally. 3)We 
detect a low ionized NLR in I ZWI, because of the low flux ratios [OIII] n /H/3 n (~ 1.7). 
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1 Introduction 



The asymmetry of [OIII]A5007 profile possessing an extended wing at blueward 
and a sharp red falloff was mentioned by Heckman et al.(1981) and by many 
subsequent researches(e.g. Busko & Steiner 1988, Grupe et al. 1999, Holt et al. 
2003, Tadhunter et al. 2001). On the basis of investigations on large sample, 
recent studies claimed that a large fraction of AGNs are of [OIII] broad com- 
ponents. Veron-Cetty et al.(2001) reported that, in general, a relative narrow 
Gaussian profile(FWHM~ 200 — 500 km s" 1 ) with a blueshifted, broad Gaus- 
sian component (FWHM~ 500 — 1800 km s _1 ) is necessary to reproduce the 
each [OIII] profile in a number of objects(24). In addition, a peak + blue wing 
structure was observed in numerous objects by Zamanov et al.(2002). The 
authors claimed that the [OIII] profiles have a narrow, unshifted [OIII] line 
with a strong blue wing. Recently, Wang et al.(2004) identified a prominent 
blueshifted, broad [OIII] component in Narrow-Line Seyfert 1 galaxy (NLS1) 
SDSS J022119.84+005628.4. 

So far, although the peak+blue wing structure of [OIII] has been detected in 
many cases, the properties of the broad [OIII] component have been rarely in- 
vestigated and poorly understood. To our knowledge, only Sergeev et al.(1997) 
analyzed the variations of [OIII] A5007 blue wing in NGC5548. The transfer 
function of the blue wing is narrow-peaked near 450 days with no response 
up to 350 days lag. One interesting question is where is the broad component 
of [OIII]A5007 from and how large is the length scale of the line emitting re- 
gion. The answers are extremely benefit to further geometrical and dynamical 
research of emission line region in AGN. 

As usual, the study of spectral variations is a powerful technology in probing 
the physics of AGN. Variability can help us to set constraints on the sizes 
of different regions of AGNs and can give information about the processes 
governing the variations. The object, IZWI, is the prototype of NLS1 and 
was discovered by Zwicky(1971). In this paper we present two new spectra 
of IZWI taken in 2001 and 2002. The variations of the [OIII] broad compo- 
nent are examined by comparing them with the spectrum taken by Boroson & 
Green(1992, hereafter BG92). They examined the emission line properties of 
a sample of 87 low redshift (z < 0.5) PG quasars. The spectra with relatively 
high S/N ratio were obtained on a number of nights in 1990 and 1991. A reso- 
lution of 6.5-7A was measured from the comparison spectra. The resolution at 
this level is adequate for identifying the evident asymmetry of the [OIII]A5007 
line. 

The paper is organized as follows. The observations and date reduction tech- 
nique are described in §2. Section 3 contains the results of spectral comparison 
and underlying implications. 
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2 Observation and data reduction 

The two spectra were taken on November 16, 2001 and December 3, 2002. 
The observations were carried out with the NAOC 2.16m telescope at the 
observatory of Xinglong and the OMR spectrograph, using a Tektronix 1024 x 
1024 CCD as detector. Two different 6001 mm -1 gratings, blazed at 5500A and 
5000A, were used in observation taken in 2001 and in 2002 respectively. Both 
grating give a dispersion of about 100 A mm -1 . The slit was oriented in the 
north-south direction for both observations and we attempted to observe as 
close to the meridian as possible. Both observations were made through a 2" 
slit which produced a resolution of ~ 6A as measured from the night sky 
lines. This set-up caused that the spectral resolution is identical with that of 
BG92 at greatest degree. In the observation taken in 2001, two exposures were 
performed. The exposure time of each frame was 2000s. The two frames were 
combined prior to extraction in order to enhance the S/N ratio and eliminate 
the contamination of cosmic-ray easily. Only one frame with exposure time 
3000s was obtained in the observation taken in 2002. 

The unprocessed frames were reduced by standard CCD procedure using IRAF 
package. The CCD reductions include bias subtraction, flatfield correction 
and cosmic-ray removal. The wavelength calibration was carried out using 
helium-neon-argon lamps taken at the beginning and end of each exposure. 
The resulting wavelength accuracy is better than lA. Two or three KPNO 
standard stars(Massey et al. 1988) were observed per night for carrying out 
the flux calibration. Each of extracted spectra was transformed to rest frame 
in terms of the redshift determined by fitting the narrow peak of the H/3 line 
to a Gaussian. 



2.1 Fell subtraction and measurements 

It is clear from the inspection of the spectra that the Fell emission is a com- 
plicating factor in measuring line parameters. In order to reliably model the 
line profiles it is necessary to appropriately model the Fell complex. The con- 
tamination of the Fell complex is well subtracted by the Fell template which 
is the Fell emission in IZWI and is described by BG92. Briefly, the template 
is a two-dimensional function of FWHM and of intensity of the Fell blends. 
The detailed procedure making the Fell template can be found in BG92. The 
best subtraction of the Fell blends is derived by searching in the parameter 
space and requires a slick continuum at blue of H/3 and between 5100A and 
5500A. From each spectrum the best-fit Fell spectrum is subtracted. These 
Fell spectra and the resulting Fell-subtracted spectra are shown in Figure 
I. The spectra displayed from top to bottom panel were observed by BG92, 
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in Nov 16, 2001 (hereafter denoted as WWHI for abbreviation) and Dec 03, 
2002(WWHII, for short), respectively. In each panel, the observed spectrum 
is displayed by the upper curve. Note that these curves are offset upward ar- 
bitrarily for visibility. The Fell flux is measured between the rest wavelength 
4434A and 4684A. The upper and lower limits of the flux of the Fell complex 
are obtained by carefully iterative experiments with a series of values of the 
Fell flux. Outside of the limits, the Fell-subtracted continuum is absolutely 
unacceptable. 



2.2 Profile modelling and decomposition of other lines 



The Fell-subtracted spectrum are used to measure the non-Fell line prop- 
erties. The first step in the modelling of the spectra is to remove the con- 
tinuum from each Fell-subtracted spectrum. The continuum is modelled by 
fitting a power law to the regions which seemed to be uncontaminated by 
emission lines. In the next step, the line profiles are modelled by multiple 
Gaussian fitting in the present work(e.g. Xu et al. 2003). The modelling is 
carried out by the SPECFIT(Kriss 1994) task in IRAF package. The fitting is 
persisted until the minimum of x 2 , the measurement of goodness of the fitting, 
is achieved. In each spectrum, the profiles are modelled as follows. The each 
of forbidden lines [OIII]AA4959,5007 is synthesized from a narrow component 
and a broad, blueshifted component (see Oke & Lauer 1979, Veron-Cetty et 
al. 2004). The atomic physical relationships, F^oqt/F^q = 3(Storey & Zeip- 
pen 2000) and A4959/A5007 = 0.9904, are used for reducing the number of free 
parameters and for improving the reliability of fitting in the modelling of both 
narrow and broad components. The narrow component of [OIII] is referred 
as [OIII] n , and the broad component as [OIII] b The H/3 profile is reproduced 
by a set of three components: a narrow peak, a classical broad component 
with FWHM~ 1000 km s" 1 and a very broad base (FWHM~ 4000 km s -1 , 
see Sulentic et al. 2000b; Veron-Cetty et al. 2001; Marziani et al. 2003), be- 
cause of the strong blue wing of the H/9 profile. The width of the narrow H/3 
is forced to be identical with the width of narrow [OIII]. The narrow Hf3 is 
referred as B.(3 Q . The contribution, including the classical broad component 
and very broad component, is referred as H/3b- Figure 2 illustrates the mod- 
elling of the three observed spectra. The label of the spectrum in each panel 
from top to bottom is the same as that in Figure 1. The observed profiles are 
represented by light solid lines, and the modelled profiles, by heavy solid lines. 
In each panel, the narrow core and broad base of each line are shown by a 
short dashed line and a long dashed line, respectively. The residuals between 
the observed and the total fitted spectrum is displayed in the lower sub-panel 
underneath each spectrum. 
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3 Results and discussion 



The three observations are compared with each other to investigate the vari- 
ability of the profile of broad component of [OIII]. The two recent observations 
span only about one year, and the observation in BG92 was performed round 
about ten years ago. Table 1 summarizes the line properties of the principal 
emission lines at all epochs. In the following analysis and discussions in this 
paper, a constant flux of narrow core of [OIII] is assumed. Column 1 lists the 
indices of the observations. The indices are same as the indices in Figure 1 and 
Figure 2. The flux ratios H/3 n /[OIII] n are listed in column 2, and the ratios 
H/3b/[OIII] n in column 3. Column 4 lists the normalized flux of broad [OIII] 
component. The velocity shifts of Gaussian peaks of broad [OIII] with respect 
to that of [OIII] narrow components are listed in column 5. A negative velocity 
corresponds to a blueshifted broad component relative to narrow component. 
The column 6 lists the flux ratio of Fell emission to narrow [OIII] , along with 
the upper and lower limits. The last column lists the calculated parameter 
RFe. RFe is defined as the flux ratio of the Fell complex to H(3. We note 
that no attempt has been made to correct for the contamination of the flux 
of narrow component of H/5, because the flux of narrow H/5 is expected to be 
constant. This definition avoids the additional errors caused by the decompo- 
sition. In addition, the WWHI and WWHII averaged parameters denoted as 
the index WWH are calculated. The mean values of the interesting parameters 
are reported in the last row of Table 1, except RFe, the ratio H/?b/[OIII] n and 
FeII/[OIII] n . Note that the FWHM are excluded from Table 1, because we 
focus attention on flux variability of the broad [OIII] component. 



3.1 Variability of broad [OIII] 

By comparing the results of BG92 and WWH, we find remarkable variations 
of profile of broad [OIII] during the period of ten years separating the observa- 
tions. The normalized flux in broad [OIII] decreased from 1.90 ±0.22 to 1.00 ± 
0.28. In addition to the flux variability, the outflow (blueshift) Aiv([OIII]) in- 
creased from —701.3 ± 44.2 km s _1 to —959.7 ± 104.6 km s" 1 . The change of 
velocity, WWH with respect to BG92, was -285.4 km s" 1 which is evidently 
larger than the errorbar estimated by the multi-component modelling. 

In Figure 3 we also show the residual profiles for all epochs as an additional 
test of line variability. The residual profiles are derived by removing the con- 
tributions of continuum, H/3 and narrow component of [OIII]. The residuals 
are also normalized by a constant flux of narrow [OIII]. The vertical lines from 
bottom to top denote the centers of the modelled broad [OIII] components. 



5 



The authenticity of the variations is verified by following two pieces of ev- 
idence. In the first instance, the normalized flux in broad [OIII] remained 
constant in the recent two epochs. Once again, the coincident is elucidated by 
the superposition of the two residual spectra in Figure 3. According to this 
coincident, we are confident that the detected variations can not be caused 
by the observational fluctuations. Secondly, the flux ratios H/3 n /[OIII] n are 
constant approximately in the all epochs, which agrees with the prediction of 
generally accepted unified model of AGN. Therefore, we infer that the identi- 
fied profile variations can not be explained by the errors caused by the profile 
modelling. 

According to the variations of broad [OIII] during about 10 years separating 
the observations, we conclude that the physical length scale of broad [OIII] 
emitting does not exceed 10 light years, corresponding to ~ 3pc. On the other 
hand, the broad [OIII] is emitted from a lager region than the BLR(> 1 lyr), 
because it has not varied during the 1 year separating the observations. Thus, 
the broad [OIII] emitting region is definitely situated at outside the BLR. 

The radius of the broad [OIII] emitting region can only be obtained by indirect, 
model dependent arguments. Assuming that the gravitational force of the 
central black hole dominates the kinematics of the emitting region of broad 
[OIII], the broader profiles, with respect to the narrow profiles, indicate that 
the material emitting broad [OIII] should be located in NLR inner regions. 
The distance from central engine to the line region emitting broad [OIII] can 
be expressed as 

i?(lt-days) = 6.7 x lO" 6 ( ^ WHM - ) ( **™) • (1) 
V 3 ' V10 3 km s' 1 ) \M Q J K J 

where Vfwhm is the FWHM of the emission profile of line emitting gas. The 
inferred black hole mass of IZWI is log(M/M ) = 7.26 ± 0.11. The value was 
estimated based upon single epoch optical spectroscopy(Vestergaard 2002). 
Adopting the averaged FWHM of the broad [OIII] component ~ 1500 km s -1 , 
the estimated distance is approximately about 100 lt-days. 

The inferred distance means that the broad [OIII] in IZWI is most likely 
to be emitted in a transient emission line region(TLR) whose distance from 
central source is of order 100 lt-days. In fact, it is a long time to realize 
that there appears to be some interconnection between the classical BLR 
and NLR(Osterbrock & Mathews 1986, Sulentic et al. 2000a). Following this 
concept, a few authors investigated the properties of the TLR. For instance, 
by examining the profile of high order Blamer series and [OIII], the TLR, 
with a scale size of order lpc and an intermediate FWHM(~ 1000 km s _1 ), 
was reported in X-ray selected AGN RE Jl 034+396 (Mason et al. 1996). In the 
same object, by reproducing the observed spectra in terms of photoionization 
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calculations, Puchnarewicz et al.(1995) argued that the density of TLR(~ 
10 7 ' 5 cm" 3 ) is much lower than that of BLR and higher than that of NLR. 

Up to now, in addition to the viral motions in gravitational potential of galaxy 
bugle(e.g. Nelson & Whittle 1996), two other kinds of mechanism, radial flow 
from nucleus(for example, Crenshaw & Kraemer et al. 2000) and radio jet 
expansion(e.g. Bicknell et al. 1998, Axon et al. 1998; Nelson et al. 2000), are 
suggested to explain the observed NLR. Both models can also give a rea- 
sonable origin of the asymmetry of [OIII]. For the first case, the observed 
peak+blue wing structure can be interpreted as the emission from a radial 
outflowing component with a speed of a few hundreds km s _1 associated with 
a disk component described by circular rotation(e.g. Kasier et al.(2000) in 
NGC4151; Veilleux et al.(2001) in NGC2992). On the other hand, on the ba- 
sis of the jet expansion assumption, Nelson et al.(2000) modelled a set of long 
slit 2-dimensional spectra. This model reproduces the observed [OIII] profile 
in IZW1, if the far aside jet is obscured by galactic gas and the redshifted 
gas associated with near side jet lobe is optically thick(see Fig. 7 in Nelson et 
al.(2000)). Cecil et al.(2002) used this scenario to analyze the deep spectra of 
NGC 1068 taken at high spatial resolution. 



3. 2 Variation of Fell emission 

In addition to the variations of the broad [OIII], the variability of the optical 
Fell emission is identified by examining Table l(see column 6). We find a 
positive relation between the Fell emission and flux of H/3b- According to 
the results of a great deal of AGN monitoring studies, generally, the H/3 flux 
closely relates to the continuum luminosity(e.g. Peterson et al. 2002, Kaspi et 
al. 2000, Kollatschny et al. 2000). Therefore, it is expected that the emission 
of Fell increases with incident ionizing continuum in I ZWI. On the contrary, a 
marginally negative relation between RFe and H/3(or continuum) flux is found 
in I ZWI. This means that in I ZWI the variations of the optical Fell blends 
are weaker in comparison to H/3 line. 



3.3 Low ionized NLR in I ZWI 

The line ratios of the NLR in NLSls are always interesting parameters. The 
study of the NLR in NLS1 is not straightforward, however. Deblending the 
optical permitted lines in NLS1 is difficult because no transition between the 
narrow and broad components is observed. There is always large uncertainties 
in determining the fraction of H/3 that is emitted by the NLR. In our analysis, 
both BG92 and WWH spectra indicate that the ratio [OIII] n /H/3 n is approxi- 



7 



mate 1.7. Such low value implies that the ionization stage of NLR in I ZWI is 
much lower than the generally adopted value([OIII] n /H/3 n ~ 5 — 10). 

Is that an universal phenomenon in NLS1? As a matter of fact, Veron-Cetty et 
al.(2004) found a very low ionized NLR in I ZWI from line system N3 in which 
most of oxygen might be in the form of non-ionization, both because of the 
weakness and absence of [OIII]A5007, [SII]AA6716,6731 as well as [OII]A7324 
and because of the strength of [CaII]AA7291,7324 lines. Furthermore, the low 
ionization stage in NLSls was supported by the studies of Rodriguez- Ardila 
et al.(2000). They found that, on average in NLSls, 50% of flux of total H/3 
is due to emission from NLR and that the ratio [OIII] n /H/3 n varies from 1 to 
5. The model analysis carried out by Contini et al(2003) demonstrated that 
such fraction of the observed H(3 line flux presents in the NLR flux in NLS1 
galaxy Ark 564. Additionally, Wang et al.(2004) gave a ratio [OIII]/H/3 n ~1.62 
in SDSS J022119.84+005628.4 which is a recently identified NLS1 from Sloan 
Digital Sky Survey. All these results imply that normal broad line Seyfert Is 
and NLSls perhaps differ in ionization level of NLR. The NLSls perhaps tend 
to commonly have a lower ionized NLR. Although we do not know whether 
it is a truth, a vast of effort, of course, will be made until the final result is 
achieved. 



4 Conclusions 



Two new spectra of NLS1 galaxy I ZWI were recently taken by us and were 
compared with the spectrum taken by BG92. The comparison allows us to 
make following conclusions: 

(1) The variations of broad [OIII] are identified during the ten years sepa- 
rating the observations. The variations indicate that in I ZWI the length 
scale of emitting of broad [OIII] does not exceed 3pc and is larger than 
0.3pc. Assuming a Keplerian motion in emission region, the material emit- 
ting broad [OIII] is possibly located in TLR, between BLR and NLR. 

(2) A positive relation between the Fell emission and rL3(or continuum) 
flux is identified in the object I ZWI. Whereases, the parameter RFe 
marginally decreases with ionizing continuum. 

(3) In I ZWI, the inferred flux ratio of narrow component of [OIII] to narrow 
component of FL3 is ~ 1.7 rather than the generally adopted value (~ 
5 — 10). Such low value suggests that there is a low ionized NLR in 
I ZWI. 
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Fig. 1. An illustration of the Fell subtraction for the three independent observations. 
From top to bottom panel, the spectra were taken by BG92, on Nov 16, 2001 (here- 
after denoted as WWHI) and on Dec 03, 2002(WWHII), respectively. In each panel, 
the top curve shows the observed spectrum which is shifted upward arbitrarily for 
visibility. The Fell-subtracted spectrum is shown by the middle curve. The bottom 
curve shows the best-fitted Fell template. 
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Fig. 2. Profile modelling for three observations. The label of each spectrum from top 
to bottom is the same as in Figure 1. The light and heavy solid line represent the 
observed and modelled profile, respectively. The narrow core of each line is displayed 
by short dashed line, the broad base by long dashed line. The residuals are shown 
by the lower sub-panel underneath each spectrum. 
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Fig. 3. Comparison of residual profiles of broad [OIII] components. The positions 
of center of modelled broad [OIII] are shown by lines from bottom to top. 
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Table 1 

The Properties of Emission Lines of I ZWI 



Indices 3, 
(1) 


H/V[OIH]n 

(2) 


H/V[OIII] n 

(3) 


[OIII] b /[OIII] n 

(4) 


At; r ([OIII]) b 
(5) 


FeII/[OIII] n 

(6) 


RFe 

(7) 


BG92 
WWHI 
WWHII 
WWH 


0.60 ± 0.32 
0.59 ±0.87 
0.64 ±0.83 
0.62 ± 0.85 


5.79 ± 0.67 
4.96 ± 1.04 
3.56 ±0.83 


1.90 ±0.22 
1.04 ±0.27 
0.95 ±0.29 
1.00 ±0.28 


-701.3 ±44.2 
-919.1 ±95.6 
-1000.2 ± 112.8 
-959.7 ± 104.6 


12.66±i;i; 
10.961?;^ 


1.98 
1.98 
2.15 









a The indices are identical with that in Figures 1 and 2. 

b The velocity shifts, in units of km s _1 , of peaks of broad components with 
respect to that of narrow components, where a negative velocity corresponds 
to a blueshifted, broad component. 
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